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1. Introduction 

In a series of seminal papers (pQ, [2] and [3]) Cachazo, Dijkgraaf, Intriligator and Vafa 
proved that a class of hyperelliptic curves of the form 

y 2 = W'{z) 2 + f{z). (1) 

where W(z) and f(z) are polynomials such that degf = degW — 2, are essential to 
study certain supersymmetric gauge theories. More concretely for the description of 
the low energy dynamics of N = 1 U{n) supersymmetric gauge theories obtained by 
deforming pure M = 2 theories by adding a tree level potential Tr for the Higgs 

field $. These curves emerge from the duality between supersymmetric gauge theories 
and string theories on deformed 3-dimensional Calabi-Yau (CY) manifolds defined by 
the equation 

W'(z) 2 + f(z)+u 2 + v 2 + w 2 = 0. (2) 

Here the polynomial f(z) is introduced to regularize the singular CY manifold 

W' 2 + u 2 + v 2 + w 2 = 

according to the process called Geometric Transition [TJ. Furthermore, the curves ([!]) 
can be generated from a special class of Seiberg-Witten curves [I] of pure Af = 2 gauge 
theories which satisfy a certain factorization property [5]. 

The curves ([I]) also appear in the study of the asymptotic eigenvalue distribution 
as n — > oo of holomorphic matrix models with partition functions 

Z n = [ e~ nTrW ^dM, (3) 



such that the eigenvalues of the nxn matrices M are constrained to lie on some path in 
the complex plane (see [6], [7] and [8]). They also arise in the study of the asymptotic 
zero distribution as n — > oo of monic orthogonal polynomials P n (z) = z n + ■ ■ -, verifying 
the conditions (see 0, [ID] and [TTj ) 

P n (z) z k e~ nW{z) dz = 0, fc = 0,...,n-l, (4) 



r 



The curves have an even number 2g(l < q < N) of branch points given by 
the odd-order roots j3 = (/?!,..., /3 2( j), of the polynomial y(z) 2 , and we may introduce 
a system of q disjoint regular curves (cuts) jj with endpoints (/^'-i, /%) which shrinks 
to a critical point aj of W(z) (W'(aj) = 0) at f(z) = 0. These cuts define the 
regularization of q double roots aj of the undeformed zero genus curve y 2 = W'(z) 2 . In 
the applications to supersymmetric gauge theories the curves are characterized by 
a set of q complex numbers s = (si, . . . , s q ) called 7 t Hooft parameters, which represent 
the gaugino condensates of glueball superfields and are given by 

f y(z)dz = -Air is j, j = l,...,q, (5) 

J Aa 
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where Aj are counterclockwise cycles around the cuts jj and y(z) is the branch of ([T]) 
with asymptotic behavior 

y(z) = W'(z) + 0(z~ 1 ), z^oo. (6) 

The aim of this work is to study the structure of the moduli space of the curves ([!]) 
where 



y N+l 



N N 



n=l k=l 

We assume that the coefficients g^ are fixed while the coefficients are considered 
as arbitrary varying complex parameters. Henceforth we will refer to ([!]) as spectral 
curves and to the coefficients t = (tx, . . . , t^) as the deformation parameters. Our main 
interest is to apply the theory of integrable systems to describe the analyticity properties 
of the roots of y(z) 2 as functions of either the 't Hooft parameters or the deformation 
parameters. 

From the point of view of the underlying U(n) gauge theory, g-cut spectral curves 
are associated to classical vacua such that the n eigenvalues of the Higgs field $ distribute 
among q critical points a, of W(z) with multiplicities rij where n\ + . . . + n q — n. On 
these vacua the gauge group is spontaneously broken into q factors 

U(n) -> U(m) x U(n 2 ) x ■ ■ ■ x U(n q ). 

The role of spectral curves in the low energy dynamics of the U(n) gauge theory 
is of central importance since they determine the prepotential functional (see PQ, [2] and 

i) 

7= I p{z)W{z)\M- \ I [ P (z)p(z')\og(z-z'f\dz\\dz'\, (8) 

which is the key object to obtain the vacuum expectation values of many important 
operators of the quantum theory. Here 7 denotes the union of all the cuts 

j = Tl u • • • U j q , 

and the function p(z) (the spectral density) is defined by 

P (z)\dz\ = ^— = -^— , zer (9) 

The subscripts ± in z± indicate the limits of functions from the left and right of the 
cuts. The 't Hooft parameters can be considered as partial charges with respect to the 
spectral density since we may rewrite ^ as 

p{z)\dz\= Sj . (10) 

n 

The prepotential is independent of the precise form of the cuts as long as they remain 
in their respective homology classes in the complex plane with all the branch points 
deleted. 

In the present work we introduce a decomposition of the moduli space of spectral 
curves onto classes in terms of sets Ai q of spectral curves with at most 2q branch points. 



Spectral curves in gauge/string dualities: integrability, singular sectors and regularization. 4 



Each class M. q is naturally stratified into a regular and a singular sector. The regular 
sectors contain the spectral curves which generically are determined by the constraints 
(J5|, so that they depend on q moduli which may be either the 't Hooft parameters s 
or a subset of q deformation parameters t. In the language of supersymmetric gauge 
theory [12], the regular sectors represent the different 

quantum phases of the vacua of the gauge model and, consequently, the singular 
sectors describe the regions of phase transitions [13J. 

One of the main properties of spectral curves Q constrained by ^ is that the 
function 



W(z) - 



[ p{z')\og{z-z'f\dz'\ (11) 

is locally constant on 7. In this work we use this property to show that the prepotential 
functional is a natural notion associated to spectral curves which do not require the 
introduction of an auxiliary random matrix problem. Thus we prove that the spectral 
density (|9| is an extremal point of the prepotential functional constrained by (10) and 



we provide a simple proof of the r-function character of the logarithm of the prepotential 
functional (J8|. 

We revisit the theory of Whitham equations for spectral curves [El [15], [16] to 
find the differential equations which characterize the classes A4 q and to determine the 
spectral curves from systems of implicit equations involving the roots of y 2 (z) and 
the moduli. Furthermore, we provide an alternative way to characterize the spectral 
curves ([I]) from of a reduced system of implicit equations for the branch points involving 
the coupling constants and the q deformation parameters (ti, . . . ,t q ) only. The reduced 
system represents the critical points of certain polynomial solutions W of Euler-Poisson- 
Darboux equations [17], and it turns out to be a generalization for multicut spectral 
curves of the planar limit equations of one-cut (q = 1) random matrix models [18J- 
PJJ]. The main advantage of the use of reduced systems is that they drastically simplify 
the analysis of singular sectors where the gradient catastrophe for the branch points 
and, consequently, the breaking of analyticity for the spectral curve and the free energy 
(|sj) happens. In this way we show how singular spectral curves corresponding to the 
emergence of higher order branch points can be described in terms of simple classes of 
degenerate critical points of the solutions W of the Euler-Poisson-Darboux equations. 
Moreover, the reduced systems leads naturally to the multiple scaling limit method 
to regularize the critical behaviour of the branch points. Thus we propose a method 
of regularization of the gradient catastrophes based on substituting the critical point 
equations for the functions W by the Euler-Lagrange equations for a density functional 
W re9 obtained by adding appropriate terms with derivatives to W. In the simplest cases 
the corresponding regularizing differential equations are the Painleve-I equation and its 
multi-component extensions. 

The paper is organised as follows. In section 2 we analyze the moduli space of the 
spectral curves and introduce the classes M. q as well as their regular and singular sectors. 
We obtain the variational characterization of the spectral density from the prepotential 
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functional. Then we discuss the integrable structure supplied by the Whitham equations 
for determining spectral curves from 't Hooft parameters. Section 3 deals with the 
analysis of the moduli space of spectral curves in terms of deformation parameters t 
and with the study of the analyticity properties of the roots of y(z, t) 2 as functions 
of the deformation parameters. We illustrate our theoretical scheme with a complete 
description of the moduli space of spectral curves for the Gaussian and cubic models 

W(z) = -, --gz. 

Section 4 presents the method of reduced systems based on solutions of Euler-Poisson- 
Darboux equations to determine spectral curves. We also show the connection between 
singular spectral curves and degenerate solutions of reduced systems. Sections 5 and 
6 describe a simple implementation of the multiple scaling limit to regularize gradient 
catastrophes of the branch point functions /3(t). Moreover, the natural emergence of 
the Painleve-I equation is demonstrated. 



2. The moduli space of spectral curves 

2.1. Regular and singular sectors of spectral curves 

We may distinguish different classes in the set of spectral curves Q by looking at the 
multiplicities of the roots of y(z) 2 

K,mi 2) ---]) (1 < zi < z 2 < ■ ■ ■) 

where mj is the number of roots with multiplicity equal to i. 

Since the polynomial y(z) 2 is of order 2N it has an even number 2q of roots f3{ with 
odd multiplicities 2q^ + 1 (branch points of y{z)) and a certain number p of roots ai 
with even multiplicities 2pi so that 

P 2g 

y{zf = \[{z - atf^ \[{z-^ + \ Pl > 1, Qi > 0, (12) 
i=i i=i 

where 

P 2q 
1=1 i=l 

If we equate ([T| and (12) and identify coefficients of powers of z we get a system of 2 N 
equations for the p + 2q + N variables (ai,(3i,tk). Thus, fixing the root multiplicities 
and assuming that all the roots («;, Pi) are distinct we obtain a family of spectral curves 
which depends generically on p + 2 q — iV parameters. 

We are interested in spectral curves determined by q constraints of the form ^ 
where 2q is the number of odd order roots. Thus we focus on families of spectral curves 
which generically depend on q free parameters or, equivalently, such that p + q = N. 
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This means that y(z) 2 has root multiplicites given by [(2q)i,p2\, so that it factorizes in 
the form 

y\z) = \{{z- ai ) 2 \{{z-l3 i ), p + q = N. (14) 

1=1 i=l 

In order to deal with the structures of singularities we will consider degenerations of the 



factorization (14) in which some roots coalesce. This motivates the following definition: 



Definition 1 Given q = 0, 1, . . . , N we denote by M q the set of all the spectral curves 
which have 2q branch points (q cuts) at most. 

The class M q represents the set of all spectral curves which admit a factorization 



of the form (14) where some of the roots (a/,/3j) may coincide. Obviously we have 



M C Mi ■ ■ ■ C M N . (15) 

Definition 2 We define the regular sector iegM q as the subset of the spectral curves 
in M q which have all the roots (a/, distinct. We define the singular sector of M q by 
singA^ = M q \ regM q . 

The regular sector regM q represents the set of spectral curves such that [(2g) 1? p 2 ] 



or, equivalently, which factorizes in the form (14) with distinct roots (aj, The 
spectral curves in the singular sector singA^ admit a factorization of the form (14) but 



such that some of the roots (ai,/3i) coincide. It is clear that M q > C singA^ g for q' < q. 
2.2. The prepotential functional 

Next we will discuss an important property of the prepotential functional in the theory 
of spectral curves. To this end let us consider a spectral curve in the regular sector 
regA'lq and take a semi-infinite oriented regular path V (see figure 1) containing all the 
cuts 7j such that for each z' in V there exists a branch of log(,2 — z') that depends 
analytically of z in C minus the semi-infinite arc T z > of F ending at z', and verifies 

\og(z + —z')+\og(z^—z') = \og(z' + —z)+\og(z'__—z), for all z ^ z' in T.(16) 



The property (16) is essential for the consistence of the following definition of 
\og(z — z') 2 

\og{z - z') 2 = \og(z + - z') + log(^_ - z'), z, z! e T, (17) 



which is assumed in the expression of the prepotential (|8j). It is easy to prove that 
for arcs F with parametrization z = z{r) such that one at least of the functions x(t) 



and y(r) is an strictly monotone function of r, the property ( 16 ) is satisfied by the 
logarithmic branches defined by 

r a 

\og{z - z') = log \z -z'\+ / z'eF, zeC\ IV, (18) 

. r . u - z' 
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Figure 1. The semi- infinite path T ending at faq- 



where T ZjZ i is any path in C \ IV connecting z' + \z — z'\ to z. For example if T is a real 
interval of the form (— 00, x Q ] then (18) determines the principal branch of log(z — z') 
and (17) implies \og(z — z') 2 = 2 log \z — z'\. 

If we assume (IT6j) then the prepotential can be written as 



T 



p(z)W(z)\dz\ 

-\ J J p{z)p{z')(\og{z + -z')+\og{z--z'))\dz\\dz'\ 



Let us introduce the function 



g(z)= I \og{z-z')p{z')\dz'\. 

J T 



(19) 



(20) 



Note that both y(z) — W'(z) and —2g'(z) are analytic in C \ 7, vanish as z — > 00 and 
according to ([9]), have the same jump on 7. Therefore we obtain the identity 

p(z> 



\dz'\ 



z — z 



(21) 



y (z) = W\z) - 2g'(z) = W\z) 

which determines the spectral curve from the complex spectral density. Moreover, using 
y(z+)+y(z-) = 1 ze 1} (22) 

we get 

W'(z)-(g'(z + )+g'(z_)) = 0, ze 1 . (23) 

This equation means that W(z) — (g(z + ) +g(zJ)) is constant on each connected piece of 
7 or, equivalently, that there are (not necessarily equal) complex numbers U such that 

W(z)-(g(z+) + g(z-)) = l i , z e 7*. (24) 

A straightforward calculation using (16), (17) and (19), taking into account ([9]), shows 
lent to the 

■F+X)M- S '- / >|<1~| ) =0. (25) 



that (24) is equivalent to the the variational equation 

1 / /■ N 



5 

Sp 
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In this way the spectral density p(z) of a spectral curve is characterized as an extremal 



point of the prepotential functional T constrained by (10). 



From the property (24) it also follows that 

i 



T=\\ P {z)W{z)\dz\ + \Y jSj l j . (26) 
1 h 1 j=1 



2.3. Integrable structure: Whitham equations 



Given a spectral curve (14) in the regular sector regAi q we may consider the genus q — 1 
hyperelliptic Riemann surface M associated to the reduced curve 

2q 

I/O 3 = 110* "A)- (27) 

i=l 

The reduced curve encodes the complete structure of the spectral curve since using the 
identity 

p 



n<*-°«>. < 28 > 

1=1 



yo{z) 



where © stands for the sums of nonnegative powers of Laurent series in z, we may 
determine the double roots ai as functions of the simple roots Therefore, in order to 
construct g-cut spectral curves obeying ^ we must prove that the branch points (3 can 
be determined as functions /3(s) of the 't Hooft parameters. We are going to use the 
theory of Abelian differentials in Riemann surfaces (see Appendix A) to find a system 
of implicit equations determining the functions /3(s), which arises as a consequence of 
the integrability of a system of Whitham differential equations. 

Let {A{, BiYiZl be the homology basis of cycles in M showed in figure 2, where 
Ai are counterclockwise cycles around the cuts •y i with endpoints (/#2i-i> p2i)- We will 
denote the corresponding periods of a differential du in M by 

Ai(du) = I du, Bi(du) = I du. (29) 



Let us introduce the meromorphic differential 

d^= l -{y{z) + W\z))dz (30) 

where y(z) denotes the function on the Riemann surface M determined by two branches 
of given by 

y 2 (z) = -yi(z) =y(z), (31) 

where 

y ( z ) = W\z) - 2 g'(z) = W\z) - — + 0(z~ 2 ), as z -> oo. 
Here s denotes the so called total t 'Hooft parameter 

S = Si H h S q . 
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Figure 2. The homology basis. 



We are going to decompose dS in terms of the set Abelian differentials of first kind 
(i.e., holomorphic) {dy?;}^ 1 , second kind dVt n (n > 1) and third kind dQo described in 
Appendix A. These Abelian differentials can be written as 



d^(z) = ^44dz, da 



n 



for appropriate polynomials Pj(z) and P n {z) 
We first observe that 



2 Vo(z) 



dz, dQr 



yo(z) 



dz, (32) 



y(z)dz = < 



W'(z) 



2s 



+ 0(z ) I dz, as z — > ooi, 



(-W'(z) + j + 0{z~ 2 )\ dz, 



(33) 



as z — > oo 2 . 



Since the only poles of y(z)dz are at ooi and oo 2 , from (33) we deduce that 

N+l 

d§ - 9n d^n + s dQ (34) 

n=l 

is a first-kind Abelian differential in M. Consequently it admits a decomposition in the 
canonical basis 

N+l q-1 

d§ — g n df2 n + s df2o = Xj d^-, (35) 

n=l j=l 

for some complex coefficients Xj G C. Thus, we may write 

q-l N+l 

y{z)dz = -W'{z)dz + 2 ^ Xj d<pj + 2^2g n dQ n - 2s dQ - (36) 

j=l n=l 

Moreover, since Ai(dVl n ) = for all i = 1, . . . , q — 1 and n > then from ^ and (36) 
we get 



A, 



-2tt i Sj. 



Therefore we obtain the decomposition 

N+l q 



d§ = ^ g n d£l n - Y s j ( dfi o + 2 7r i (1 - S jq ) d<pj) 

n=l j=l 



(37) 
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The differential dS generates a system of Whitham equations with respect to the 't 
Hooft parameters s [15]- [16]. To derive this system we notice that d§ has poles at ooi 
and 002 only. Moreover its derivatives with respect to the 't Hooft parameters satisfy 

' < 1 



ddS 

dsi 



+ 0(z ) dz, as z — > ooi, 



+ 0{z 



-2 



and 



which imply 



A, 



dd§ 

8Sn 



-2 7ri (5 



■jk, 



dz, as z —>• oo 2 . 



q-1, 



ddS 

dsi 



-dfl — 2 7r i (1 — 5j g ) difj. 



As a consequence we have the following system of Whitham equations 



d_ 
dsi 



(dn + 2 7T i (1 - <y d^J = — (dn + 2 7T i (1 - s iq ) d<p t j . 



Using (|32j) we may rewrite (|38|) as 

dy{z) 



P Q (z)-27ri(l-5 jq ) Pj (z) 



(3f 



(39) 



(40) 



dsj y (z) 

Hence dividing both members of the resulting equations by yo(z) and equating the 
residues at z = Pi we get the following system of ordinary differential equations 



<9A 
d Sj 



P o (P i ,0)-2iri(l-6 jq )p j (p i >0) 



(41) 



nr=i(A - a i) n^i(A - fa) ' 

which will be henceforth referred to as the Whitham equations for the branch points. 
From these equations it follows that for spectral curves in the regular sectors the branch 
points are analytic functions of the 't Hooft parameters. 
From (36) and using (32), we have 

N+l q-1 

y(z)y (z) =2j29nPn(z)-2sP (z)-4niJ2s j p j (z). (42) 

n=l j=l 

Therefore if we set z = Pi in this identity we find the system 

N+l q 

Y,9nPn{Pi,P) -Y i s j (p (p i ,f3) + 2iri{l - 5 i? ) Pi (A,/3)) = 0, (43) 

n=l j=l 

where i = 1, . . . , 2q, which provides an implicit integration of the Whitham equations 
and determines the branch points in terms of the 't Hooft parameters. We observe that 
(43) means that the generating differential dS vanishes at the branch points z = Pi 
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2.4- Prepotentials and T-junctions 



As we have said before, the present work deals with the moduli space of spectral curves 
((!]) for a fixed polynomial W(z), but obviously, the coefficients g = (gi, . . . , g^) of W(z) 
can be also considered as free complex parameters too. In that case, repeating the 
analysis of the last subsection we immediately find 
<9d§ 



dg n 



dfL 



and consequently we may enlarge the systems of Whitham equations (39) and (41) 
including differential equations with respect to the parameters g. In this sense it is 
important to notice the properties of the prepotential functional as a function of these 



parameters. Indeed, from (|8|), (19) and (24), taking into account that p(/3j) = 0, we 
have 

()T ' z n p(z)\dz\ + I [W(z) 



dg n 



Jlogiz-zTpiz'W^^ldzl 



z n p(z)\dz\ + / (W(z)-(g(z + )+g(zj)) 



dpjz) 
dg n 



\dz\ 



I z n p{z)\dz\+j2h I 
J l j=l J Tj 



3=1 

q 



dpjz) 
dg n 



r H Ft* 



|dz| 



z n p(z)\dz\. 



On the other hand from ( |21[ ) and (|30j) it follows that as z — > oo\ 

d§ = (w'(z) - E ^ / V)^W)M)te. 



n>l 



Then we conclude that 

dT 1 



z n d§, 



(44) 



dg„. 27r i 

where is a positively oriented circle of large radius. 

It is immediate to see that repeating the same procedure for the derivatives of the 
prepotential with respect to the 't Hooft parameters one obtains 
cXF 

dsi 



(45) 



Moreover, using (24) and (31) we find that the Lagrange multiplers li are related to the 
5-periods 



I dS=(l q -h), 

JBi 



so that using ( 45 ) we have 
ds n dsi 



dS, 



,q- 1, 



l,...,q 



(46) 



(47) 



Bi 



The relations (44) and (47) mean that the logarithm log J 7 of the prepotential functional 
is a r-function of the Whitham hierarchy with respect to the variables (g, s). 
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3. Spectral curves in terms of deformation parameters 



In the previous section the set s of 't Hooft parameters has been used to describe the 
integrable structure which underlies spectral curves. However, the use of the alternative 
set t of deformation parameters turns to be more convenient for dealing with several 
problems as, for example, the construction of explicit examples of multicut spectral 
curves. 

3.1. Roots of spectral curves as functions of deformation parameters 
The spectral curves in A4 q satisfy the system of 2N equations 
P(a l ,t) = P'(a l ,t) = 0, l = l,...,p, 



(4£ 



P(A,t) = o, 

where P(z,t) denotes the polynomial 



l,...,2q, 



N 



N-k 



(49) 



P(z,t) = y(z,t) 2 = W\zf + J2tk 

fc=i 

This system defines an affine variety in the space 

C 2N+q of variableg ( ajj The 

ie system (48) 

d(p(ai), P(a p ), P'{ ai ), P'(a p ), P(f3 1 ), • • • , P{fr 



Jacobian of the system (48) 



2g. 



d\a,i, ■ ■ ■ , a p , /3 2g , ti, ■ ■ ■ , t^j 



(50) 



evaluated at solutions of (48) is of the form 



A 


B 


C 


D 



(51) 



where A is the p x {p + 2q) zero matrix since from (48) 
Ai m = 5i m P'(ai) = 0, l = l,...,p, m = 
Furthermore B is the Vandermonde matrix 



,p + 2q. 



( < 



N-l „N-2 



B 



a 



JV-l „N-2 

Otr, 



1 \ 

1 

1 J 



C is the diagonal matrix 

C = diag(p"(a 1 ),---,P"K), P\^) 
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and 



/ (iV-l)af- 2 (N-2)a{ 
(N-l)a 



D 



N-3 



N-l 



(N - 1)< 



Pi 

p 



N-l 



(N - 2)a»- 



(N - 2)a»- 



N-l 



Pi 

P 



N-2 



N-2 



iN- 
>2q 



o\ 




1 
1 



/ 



V #T Pi 

Let us consider now spectral curves in the regular sector regA^ g . It is easy to see 
that for any set of q deformation parameters (i^, • • • , tj ) the minor of the Jacobian 
corresponding to the derivatives with respect to the variables (otj, Pi) and the remaining 
N — q deformation parameters tk (k ^ {ii, . . . , i q }) is a non singular (2N) x (2N) matrix. 
Therefore the implicit function Theorem implies 

Theorem 1 For spectral curves in the regular sector regA4 q , each set of q different 
deformation parameters (tj i; • • • , t iq ) defines an analytic system of coordinates, i.e. all 
the variables (a>i, tk) are local analytic functions of (t^, ■ ■ ■ , t iq ) 

If we eliminate the variables in the system (48) or, alternatively, in the 



system obtained identifying coefficients of powers of z between ([T]) and (14), we get a 
set of N — q constraints 

/*(t) = 0, k = l,...,N-q, (52) 

which characterizes the class Ai q in the space of deformation parameters. These 
constraints can be used to define analytic systems of coordinates in regA^g given by sets 
{Ud • • • j Uq) °f Q distinct deformation parameters. 

The class Mn contains all the spectral curves ([TJ and will be called the generic 
class. Its elements can be factorized in the generic form 

2N 



2/0, t) 



(53) 



where the roots /?, may coincide. In the case of the generic class the set of the 
deformation parameters (ti, . . . ,t^) defines an analytic system of coordinates in the 
regular sector regA^Ar. Obviously the domain in the space of deformation parameters 
which corresponds to reg.M n is the open set determined by 

D(t) ^ 0, (54) 

where D = D(t) is the discriminant of the polynomial y(z,t) 2 . Consequently singA^jv 
is represented by the hypersurface determined by the constraint 

D(t) = 0. (55) 

From the point of view of the applications to gauge/string dualities it is important 
to determine under what conditions a map (t^, • • • , t iq ) (si,...,s q ) relating the 
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coordinate charts of deformation parameters and 't Hooft parameters defines an analytic 
change of variables. Now from (|5| we conclude that the Jacobian matrix of these 
transformation is given by 



dt ik 87ri/ A A rf! A ^., dt ik Jy{z) 



Thus the question reduces to the verification of the nonsingularity of this Jacobian 
matrix. For spectral curves in regA'f n it follows from ^ that 

1 dz, j,k = 2...,N, 



dtk 8ni T A . Ftt2N 



This means that the Jacobian of the transformation (t 2 , • • • , t/v) ( s 2, • • • , sjv) is non 
singular since it is the matrix of A-periods of a basis of Abelian differentials in the 
Riemann surface given by 

2N 

y(*) a = ![(*-&)> 



where all the branch points are different (N — 1 genus). Moreover from (21) and (10) it 
is clear that 

t x = -As. (56) 

Therefore the change of variables (ti,...,t N ) \-¥ (s%,...,sn) is an analytic 
diffeomorphism locally around any point in the regular sector reg.Mg. Consequently 
using Theorem 1 we conclude that the set of 't Hooft parameters (si, . . . , Sjv) provides 
an analytic system of coordinates in regA^iv. 

3.2. Whitham equations with respect to deformation parameters 

We have seen that for spectral curves in the regular sector regA^ g all the variables 
(ai,Pi,tk) can be expressed as analytic functions in any coordinate system of q 
deformation parameters. In particular one may consider the expressions for the simple 
roots in terms of the coordinates . . . , t q ) 

A = A(ti,t 2 ,...,g. (57) 
We are going to derive the system of Whitham equations which characterizes the 



functions (57). Indeed, taking into account that the double roots on are distinct then 
from (49) we get the identities 

p i 
ti+« = Yl V ^ a ) [ p ( a m,t) ~ (W\a m f + Y,tj , (58) 

m=l j=l 

for I = l,...,p = N — q where Vi m (a.) is the inverse of the Vandermonde matrix 
(af~ m )f m=1 . As a consequence the functions 

p 

p(a, t) - J2 v ^ a ) p («'- *) ( 59 ) 

l,m=l 
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do not depend on the deformation parameters . . . , tjv) since they can be written 



as 

p 



Q(a, ft) + J2 *i ~ E /T 1 M«) 



(X 



j"=l (,m=l 

where 



l«™) 2 . 



i,m=l 

Now for spectral curves in regA4 g we may express all the variables (az, tfc) as functions 
of the coordinates (ti, . . . ,t q ). Hence differentiating the equations 

P(A, t) - E ^ ^«(a) P ( a m, t) = 



with respect to the free parameters {i,}J = i, taking into account that P(a/,t 
P'(a;/,t) = 0, we obtain the system of differential equations 

dh z kn i*r l vu<*)<- i -tf- i 



(60) 



This is the system of Whitham equations for the branch points with respect to the 
deformation parameters. 

The system (60) shows how the functions Pi(t\, . . . , t q ) become singular for spectral 
curves in the singular sector sing.Mg, where at least two roots of the set (a;, /%) coalesce. 



For q = N and q = N — 1 the system reduces to 
djk _ 

% n^(A-/3i)' 

and 



(61) 



^ + - / flVJTT / O ftV x a=-Ng N (62) 

respectively. 

5.5. T/ie Gaussian and cubic models 

A complete explicit analysis of the classes A4 q is possible only for polynomials W(z) of 
low degree. We next discuss some examples for the cases of degrees two and three. 



Gaussian model For the Gaussian model 

W(z) = £ (63) 



we have N = 1 and y(z,t) 2 = z 2 — 4s. The system (48) for q — 1 and g = reduce to 
P 2 -4s = 0, z = l,2, 
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and 



ct 2 -4s = 0, a = 0, 



respectively. Thus for Ai\ the branchpoints are 

= = 2^- (64) 

The regular sector regA^i is given by the spectral curves corresponding to values s ^ 0. 
The singular sector singA^i coincides with A^o and contains only the spectral curve 
y(z) 2 = z 2 which corresponds to s = 0. 

Cubic model Let us consider the cubic model 

W(z) = j-gz, 9 GC\{0}. (65) 
In this case N = 2 and 

y(z,t) 2 = (z 2 -g) 2 + t lZ + t 2 , t x = -As. (66) 
The generic class M 2 is the set of all the spectral curves (J66j) . If we write 

y(z, t) 2 = {z- fa)(z - 2 ){z - (3 3 )(z - ft), (67) 

introduce the variables 

ft +ft /3 3 + /3 4 

(68) 

, , fa-fa 

°1 = o > °2 = o • 



and identify the coefficients of powers of z in equations (66) and (67) we obtain the 
relations 



u 2 = -u 1 , 5l = g -u\ + — , 5 2 2 =g-u\-—, (69) 

and the equation which determines u\ as a function of the deformation parameters (s, t 2 ) 

Au\ - Agu\ - t 2 u\ - s 2 = 0. (70) 

The regular sector regA^2 is determined by the points (s,t 2 ) G C 2 at which the 
discriminant of y(z, t) 2 does not vanish, i.e. 

regM 2 o {(s, t 2 ) e C 2 : 27s 4 + (18gt 2 + 16c/ 3 ) s 2 - t\ - g 2 t\ ^ 0}. (71) 



The singular sector singA^ is given by the spectral curves (66) which have a multiple 
root. Their possible sets of root multiplicities are 

[2 1; 1 2 ], [li.lg], [2 2 ], [1 4 ], 
which obviously means that 
singA^ 2 = Mi. 
Let us now consider the class Mi 

y ( z ,t) 2 = (z-a) 2 (z-fa)(z-fa). (72) 



Spectral curves in gauge/string dualities: integrability, singular sectors and regularization. 17 



Equating (66) and (72), identifying coefficients and using the variables 

Pi +fa , 01 ~ 02 

we get the system of equations 

a = -u, 2u 2 + 5 2 = 2g, 
u5 2 = 2s, u A -u 2 5 2 -g 2 = t 2 . 1 ' 

Therefore u satisfies the cubic equation 

u 3 - gu + s = 0, (74) 

and S 2 is determined by 

5 2 = —, ifu^O, (75) 
u 

5 2 = 2g, if u = 0. (76) 



Moreover, eliminating u and 5 in (73) we obtain a constraint for the deformation 
parameters (s,t 2 ) which, as it should be expected, is equivalent to the vanishing of the 
discriminant of y(z,t) 2 

MiH{(s, * 2 ) e C 2 : 27s 4 + (180*2 + 16# 3 ) s 2 - t\ - g 2 t\ = 0}. (77) 



According to (73) 0\ = 02 if and only if (s, * 2 ) = (0, 0), i.e. y 2 = {z 2 — g) 2 . Moreover 
if a coincides with 0i or 2 then the polynomial y(z) 2 acquires a root of multiplicity 
m > 3 and therefore it satisfies [20] 



4 

5](-i)W(^ 2 4 -yw = o, 

k=0 



which gives the constraint 

3t 2 + 4g 2 = 0, (7i 



or, due to (77) equivalently, 



(--/-,)= ( ±2('^ 3/2 ,4<A (79) 



.3/ 3 

Therefore, the singular sector singA^i is given by only three points of the space of 
deformation parameters 

singM 1 o{( S ,* 2 ) = (0,0), (±2(J) 3/2 ,~/)}. (80) 

The corresponding root multiplicities are [2 2 ] for (0, 0) and [li, I3] for (±2 (g/3) 3//2 , — 4g 2 /3). 
Figure 4 exhibits an example of a curve in singA^i in which a triple root is formed as 
the outcome of a merging of a simple root and the double root. 

From (74) it follows that the class M.\ decomposes into three subsets 

M 1 = U 2 k=0 M[ k) 
which are determined by the three branches 

u k (s) = -^j--A k (s), (A; = 0,1,2) (81) 
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- 



-1 - 




-2 - 



Figure 3. The classes of spectral curves of the cubic model for g — 1 in the space 
of real values of (t 2 ,s). The whole real plane describes M.^ > the curve corresponds 
to the implicit equation (77) and represents singA^2 = -Mi, the three marked points 
represent singA^i. The point (0,0) represents AAq. 



where 



The expansions of Uk(s) as s — > are 

Ms) = -V9~^ + 0j2+ °( s3 )' (») 

^) = v^-|^ + ^ 3 ), (84) 

u 2 (s) = - + 0(s 3 ), (85) 
9 

which show that in this limit the branches Uq and U\ represent families of one-cut 
spectral curves in which the cuts shrink to the corresponding critical point —y/g and 
+y/g of W(z). The solution u 2 (s) has a different behavior as s — > 0, because in this 



limit 5 — » 2g, i.e., the solution reduces to (76). 



Finally, the class A4q contains only one spectral curve y 2 = (z 2 — g) 2 which 
corresponds to (s,t 2 ) = (0,0). 

4. The Euler-Poisson-Darboux equation and the reduced systems 



As we have seen the regular sectors regA^g are completely characterized by the functions 
(57). We will now prove that these functions can be determined from a system of 2q 
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Figure 4. The functions afo) and ft ^2) at a merging of the double root and a simple 
root in the cubic model for g = 1. 



equations for the variables 

These systems will be called reduced systems. They are determined from a particular 
polynomial solution of the Euler-Poisson-Darboux equation E(—-^, — -A (see e.g. |17j ) 



2 Wj - A) 



5 2 



(86) 



4-1. Reduced systems 
Let us define the function 



W(/3,t) 



1/0(2, /3) 2/(2, t) 



dz 
2vri' 



37) 



where 



2g 



JJ^-A), y(2,t) = 
\ i=i 



2/o(^, /3) 

and Too is a positively oriented circle of large radius. It obviously satisfies the Euler- 
Poisson-Darboux equation (86). Moreover, using the expansion of y(z,t) at z — > 00, 



i.e. 



t) = vWTTM) = w'{z) + + 0^ 



we deduce that W reduces to 



W(/3,t) 



1 



W'(z)y (z,(3) 



dz 
2ttT 



(88) 



W'(z) 2vri 2 ' 
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Observe that W — t q+ x/2 does not depend on . . . , t/v) and that it is linear in 

(ii, . . . , t q ). 



Examples: In the cubic model we have for q = 1 
1 



and for q = 2 



128 
1 

256 



- (fa - fa) 2 (5/3i + Qfafa + 5/3 2 2 - IQg) - 32 (fa + &) h + 64t 2 



2 (A + & + & + fa) (fa - /3 2 ) 2 (& - &) 2 



+ (A - fa) \ - 7/3? - 9fapf + hfaft + hfaft - 9(3$ fa 
+Qfafafa + Qfafafa - 7/3 2 3 + h(3 2 2 fa + h(3 2 2 fa 
+ (fa -fa)H- 7/3f + 5/3i/3f + 5/3 2 /5f - 9/3 4 /3f - 9/3 4 2 /3 3 



+6/W3 + 6fafafa - 7 fa + 5/?!/?? + bfafa 



+ T7:(fa + fa-fa- fa) (fa -fa + fa- fa) (fa -fa-fa + fa) 9 
lb 



+ 



(fa -fa + fa- fa) 2 + 4 (fa fa + fa fa) + 8g 



(fa + fa + fa + fa) t 2 



We are now ready to prove the following characterization of the functions (57) 

Theorem 2 For spectral curves in Ai q the functions (3 = f3(t^) satisfy 
dW 



dfa 

Proof One has 

dyo(z,(3) 



0at/3 = /3(t^), i 
iyo(z,P) 



l,...,2q. 



19) 



so 



dfa 
<9W 



2 z - fa ' 
Vo(z, (3)y(z,t) dz 



dfa 2 J r ^ z-fa 2vri 



Moreover for spectral curves in A4 q we have that 

y (z,(3)y(z,t) -A- 
—77 = 11(* - «0 11^ - &)> 



* - A 

and then the statement follows. 



(90) 



The equations (89) can be expressed as 

J roo y (z,(3) V V{) + 2 W'(z) 



dz 

2~7Ti 



l,...2q, 
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and therefore for distinct Pi they are equivalent to the system 

K ' 2 W'(z) ) 2ni 



k k = 0,...2q-l 



(91) 



We will refer to (89) as the reduced system for the 3g variables (/3,t^ 9 ^). It must be 



observed that the previous system (48) for (cfy, involves a larger number 2p + 3g 
of variables . 



Examples: For q — 1 the reduced system (89 ) is the well-known pair of equations which 



determines the planar limit of random matrix models 

r w(z. 



W\z) 



dz = 0, 



(92) 



v /(z-/3 1 )0z-/3 2 ) 



dz = 47ris. 



For g = 2 the reduced system (89) takes the form 

W'(z) 

K J dz = 0, 



W'{z) 



-.dz = 0, 



W'(z) 



(93) 



ITU* -to) 



--dz = 47ris, 



W'(z) 



1 4 

dz = 4 7T i /5j - iV ) s - vr i t 2 . 



4-2. The classes TZ q of solutions and their singular sectors 

The reduced systems allow us to establish a simple characterization of relevant subsets 
of singular spectral curves. 



Definition 3 We define lZ q as the set of all the solutions (/3, t^) of the reduced syst 



em 



(89) such that fa ^ f3j for all % ^ j 
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Obviously for spectral curves in the regular sector regA^ g we have that (/3(t^), t^) G 
lZ q . Nevertheless, we are interested in degenerate solutions ((3, t^) G lZ q at which the 
Jacobian of the system (89) vanishes. 

d 2 W(/3,t) 



0. 



(94) 



From the Euler-Poisson-Darboux equations (86) it follows that if (/3, t^) e then 
all the derivatives of W at (/3, t^) can be expressed in terms of linear combinations of 
derivatives of the form c^.W. In particular if (/3,t^) e 7£ 9 we have 

d 2 W(/3,t)i ^-d 2 W(/3,t) 



n 



Thus, it is natural to define singular sectors in lZ q as follows 

Definition 4 Given a set of 2q non negative integers n = (ni, ri2, ■ ■ ■ , n2q) with some 
Hi > 1, the singular sector sing7?.g(n) is defined as the set of points (f3,t^) e lZ q such 
that 

f d k W 



fori = 1,2, 



dp? 

d n ' +2 W 
2q. 



0, VI < k < rii + 1, 



(95) 



Theorem 3 If a spectral curve in sing.Mg has roots Pi with multiplicities 2ni + 1, then 
(/3(tW),t^) e sing^(n) . 



Proof From (|87|) and (|90|) it is clear that at (3 = (3(t^) each derivative <9^ +1 W is 

IlLi(*-«0IW*-&) d ^ 



proportional to the integra 

yo(z,P)y(z,t) dz 



>-A) fe+1 2vri 



2ni 



(96) 



v 1=1 j& 

Then the statement follows since, according to our assumption, the function njLiC 2- a i) 
has zeros of order at z = Pi. 

As a consequence of this theorem we conclude that the sets sing7^ (? (n) represent 
singular spectral curves in which double zeros a\ of y(z) 2 coalesce with simple zeros Pi. 
Figure 4 exhibits a process of this type in the cubic model. 
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5. Critical behaviour and the Multiple scaling limit method 



The observation made in the previous section drastically simplify the analysis of the 
singular behaviour of spectral curves near branch points of multiplicity higher than one, 
as well as its possible regularization. 

Let us consider a spectral curve such that at a certain value t ^ the roots /3jo have 
multiplicities 2n» + 1, then at ((3(t^) } t^) we have 

d k W «9" i+2 W , s 

— tt =0, Vl< k < rii + 1, to- ^ 97 

df3 k > - - Qpn i+ 2 f V ) 



for i = l,..., 2q, and in virtue of (86 ) we also have 



0, l<k<n t + l. (98) 

To analyze the behaviour of (3(t^) near t ^ = (i 01 , . . . ,t Qq ) we will use the standard 
technique of the multiple scaling limit (see e.g. [IE]). So, we take 

tk = t ok + e a tl A = Ao + £ 7< /3*, (99) 

where e C 1 and a and 7$ are positive parameters to be determined. Now we expand 
W(/3,t (9) ) near the point ((3 ,t ( q) ). Writing W(/3,t (ff) ) = W^(/3) + ELi^O 9 )' 
taking into account (97), (98) and keeping the leading terms in the expansion, one gets 



q q 2q 

W(/3, tW) = Wo + e Q J] ^* + £ ° E E £lH t A kiPi 

fe=i fe=i i=i 

fc=l Z,m=l fc=l 



where 

W = W(/3 ,t^), A fc = V fc (/3 ), A kl = — ±(J3 0J , 

opi 



1 <9V fe 1 <9 nfc+2 W ( 3 ) 

^ m= 2^~~ (/3o) ' 5fc= K + 2)! ,9/^ +2 (/3 °' to,?) ' 



Consequently, 



df3t 

' 1 k=l k=l m=l 



+s^ +2 ^(n l + 2)B l (/3f) ni+1 + .... 



(100) 
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In the generic case when at least one A k i does not vanish, the second term is sub dominant 
with respect to the first one. In such a case the balance between the first and third terms 
is achieved if a + 7/ = (nj + 2)7;, i.e. 

7J = -^T, Z = l,...,2g. (101) 
ni + 1 

Thus, in the generic case the expansion of W(/3, t^) is of the form 

9 2<Z 



W(/3, tW) = Wo + e ^ + ^ + ^(AT 4 " 2 ) + . . . (102) 



fc=i /-i. 
where ^ = X/% =j ^^ fc j an d we put a = 1. 



Equations ( 100 )-( 101 ) imply that when and /3 are approaching to and /3 then 



9 1 

Pi-Pm~ (22(tk-t 0k )A ki y ,+ \ Z = l,...,2g. (103) 



k=\ 

and hence the derivatives of f3i blow up as 



k=l 

ast^->t^. 



m " 

dt 



1 'X)(f*-ibfcMu) (104) 





An unbounded increase of derivatives of /?; while approaching the singular sectors 
is the origin of the break of analyticity for the spectral curve and the free energy ([8]). 
From the point of view of the hydrodynamic type equations associated to the differential 



equations (60), in particular for the Burgers- Hopf equations 



wr^' 3 ^ l=1 -" N > ^'= l - Ar - l > ( 105 ) 



associated with the equations (61 ), an unbounded increase of derivatives of (3i represents 
the classical gradient catastrophe (see e.g. [2H [22])- Within the theory of polynomials 
with multiple roots and singularities of A^+i type, this singular behaviour is the 
analytical manifestation of the catastrophe arising in the transitions between two 
different strata (see e.g. [2311231125]). 

6. Regularization of gradient catastrophes 

In all physical models where such catastrophe with derivatives happens it is regularized 
at the end by one or another mechanism. It is, for instance, dissipation or dispersion for 
hydrodynamical models (see e.g. (|21j)), quantum corrections or higher genus expansion 
in matrix models (see e.g. ([18])). Formally such mechanisms consist in the appearance 
of higher order derivatives in the basic relations or differential equations of the type 



(60) or (105). Concrete forms of these higher order corrections typically are obtained 



by appropriate expansion or limit for the known model exact relations or equations. 
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Here we suggest an approach for regularization which goes in the direction opposite 
to the standard ones. The idea is to deform the scheme presented above in such model 
independent manner that it naturally leads to the appearance of derivatives in equations 
(100) and (103) to prevents the blow up of d tk fa. The principal point is the substitution 
of the equations (89) or (102) for critical points of the functions W by Euler-Lagrange 
equations 



SW reg 
5fa 



0. 



l,...,2q, 



(106) 



where W re9 is an appropriate modification of W obtained by adding terms with 
derivatives. Thus we rewrite (102), dropping the insignificant second term, as 

W(/3,t (9) ) =Wo + ey% s F I U,G8n + --- (107) 



i=i 



where Uz(/3f ) = + Bi (A*) ni+2 . Modification of the polynomial ( 107] ) is achieved by 
adding terms with derivatives of the functions fa of appropriate order in e. The general 
case is pretty involved and will be studied in a separate paper. Here we will consider 
some specific regularizations for particular situations when the corresponding equations 
(106) are rather simple. 



In the case when all n\ are different the terms in (107) are of different degrees, so 
it is natural to treat each term in the sum (107) separately. The simplest nontrivial 
modification of W(/3, t^) is of the form 

2q 



W re9 ({3*, £) =W + sY^ (e^U^f) + e Sl Q ( 

i=i 



+ 



with some constants C\ and appropriate S\. A balance of degrees in the Euler-Lagrange 

valently 

dW re9 



equation (106) or, equivalently, in the equation 

SW re9 



d t dW re9 



dfa 



o. 



implies that 



Thus, 



S, 



m + 1 



2 i 

w e9 ([3*, $ =w + J2^ + Bi (/?D ni+2 + a(%A*) 2 ) + • • • (108) 

1=1 



and the Euler-Lagrange equations take the form 



2Q 



a 2 A* 



i l + {m + 2)B l {fa l 



*\ni+l 



.,2q. 



(109) 



In the case when ni = 1 for one of the equations (109) then the simple redefinitions 



8C, 



2s I 



B, 



X 
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convert this Z-th equation into the Painleve-I equation 

d 2 n 



. )r2 60' + *. (110) 

This equation has appeared as the regularizing equation of different types of catastrophes 
(see e.g. [El EH]). In our case the Tritronquee solution of the Painleve-I equation (110) 
describes the regularization of the blow up (104) with n\ = 1, and the analyticity 
breaking associated with the formation of a third-order root from the merging of a 
double root a and the simple root fy. A similar situation takes place for n\ > 2. 

In the case when some of rii coincide, say, n\ = n 2 = ... = nk = n the polynomial 

W is 



K 

W(/3, t^) = W + + B ^iT +2 ) 

i=i 



111) 



2q 



i=k+i 



The dominant contribution of /3]J into W is of the same order. So, W re9 (f3* , £) 

should naturally contains a mixture of derivatives of • ••,/?£■ Hence, the simplest 



natural ansatz for the modified W of (111) would be 



W eff (/3*,£) = Wo + £^ (fcfl + 5/(/3D n+2 ) 
1 A <9/3* l * 88* dB*i 



m,p,q=l 
2q 



m,p,q,t=l 



e £ e^i + S,(AT !+2 + C, (^) *)+••• 



112) 



i=k+i 

where Ci, D mpq , D mpqt are some constants such that D mpq 
The Euler-Lagrange equations corresponding to the variables 81, . . . , 8* k are 



Dpmqi D m pqi C *qtmp- 



m,p,q=l 



E ^TST = 6 + (" + 2)^(A*) n+1 , 

m,p=l r 



;ii3) 



while those corresponding to the variables 8* k+1 , As are the equations (109). 

In the simplest case k = 2 and under the constraints D impq = 0, the system (113) 
is equivalent to 



83* 

J± = D 122 i + Cr, + {n + 2)B 2 8* n+ \ 
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;ii4) 



d 12 i£ - Aq - 



n 



2)Bi/3J nH 



where the variables £ and rj are defined as £1 = -Dm£ + A 77, £ 2 = D\ 22 ^ + C ?7, with A and 
C being arbitrary constants. For n = 1 and A = C = it is the special Riccati system. 
For n = 2 and A = C = it is a particular instance of the differential cubic systems 
which can be viewed as the multi-component generalizations of the Abel equation (see 
e.g. |2Z|). 

At k = 2 and under the constraints D mpq = 0, with the change of variables £i = £+77 
rj and restricted to the subspace rj 

d 2 (5{ 



6 = £ 



where 



Ah/3! 



*71+1 



A 12 /3 2 *" +1 



A 21/ 3r +i + A 22 /?r +1 + A 2 i 



- the system (113) takes the form 



A{n + 2)B 1 B 22 
B n B 22 — B 2 



A 



12 



12 



4(w + 2)£ 2 5 12 
B\\B 22 — 



4(£? 22 — gig) 
B\iB 22 — B\ 2 



115) 



.4 



21 



4(n + gjgigig 
B\\B 22 — B\ 2 



A 



22 



4(n + 2)5 2 5 



11 



B\\B 22 



R 2 



A 2 



4(5 



11 



B 



12/ 



B\\B 22 — Bf 2 ' 



5, 



Zm 



Zimj • 



=1,2 



The system (113) for arbitrary k also admits special reductions under which it takes 



a form similar to (115). The systems (113 )-( 115 ) describe regularizations of particular 
catastrophes associated with the singular sectors sing7?.(n) such that n\ — . . . — — n. 
Detailed study of the properties of such systems and their implications for the string 
and matrix theory will be presented in a separate publication. 
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Appendix A 

In this appendix we briefly discuss the elements of the theory of Abelian differentials in 
Riemann surfaces that we use in section 2. 
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Let M be the genus q — 1 Riemann surface associated to (27), it is determined by 
two sheets 

M = M 1 UM 2 , M t = {Q = (y^(z),z)}, (116) 

where 

yoA z ) = -l/o,i (*), Vo,i( z ) ~ z<? > z-^oo- 
We introduce the following Abelian differentials in M: 

(1) The canonical basis of first kind (i.e., holomorphic) Abelian differentials {dipi}?^ 
with the normalization Aj(d^) = 5ij. These differentials can be written as 

<W*) = ^dz, (117) 

where the Pj{z) are polynomials of degree not greater than q— 2 uniquely determined 
by the normalization conditions. 

(2) The second kind Abelian differentials dQ n (n > 1) whose only poles are at <x>i, 
such that 

dfi„(Q) = (nz n - 1 + 0(z~ 2 ))dz, Q ->• ooi, z = z(Q), (118) 
and normalization Ai(dQ n ) = (« = 1, . . . , q — 1). It is easy to see that 

dn n = ( V- 1 + dz, (ii9) 



.2 J/o(2 

where the P n {z) are polynomials of the form 

q-2 

P n {z) = ^(z«- 1 yo (z)) e + Y,c m z l (120) 

i=0 

and the coefficients c n i are uniquely determined by the normalization conditions. 
(3) The third kind Abelian differential dVt Q whose only poles are at ooi and oo 2 , such 
that 



dn {Q) 



-^0(z- 2 )Jdz, Q^ooi 

(121} 

1 - 



z 



— + 0(z~ 2 ) )dz, Q -> oo 2 , 



and normalization Ai(dQo) = for alH = 1, . . . , q — 1. It follows that 

dQ = ^dz, (122) 
where Pq{z) is a polynomial of the form 

q-2 

P (z) = (z-'yoiz))^ + ^ (123) 

i=0 

and the coefficients coi are uniquely determined by the normalization conditions. 

We notice that all the polynomials Pj(z) and P n (z) also depend on the set of 2q branch 
points /3, but for simplicity this dependence will not be indicated unless necessary. 
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